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We adopt the Dirac model for quasiparticles in graphene and calculate the finite temperature 
Casimir interaction between a suspended graphene layer and a parallel conducting surface. We find 
that at high temperature the Casimir interaction in such system is just one half of that for two 
ideal conductors separated by the same distance. In this limit single graphene layer behaves exactly 
as a Drude metal. In particular, the contribution of the TE mode is suppressed, while one of the 
TM mode saturates the ideal metal value. Behavior of the Casimir interaction for intermediate 
temperatures and separations accessible for an experiment is studied in some detail. We also find 
an interesting interplay between two fundamental constants of graphene physics: the fine structure 
constant and the Fermi velocity. 
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I. INTRODUCTION 



Nowadays, graphene does not require any long introduction. Its exceptional properties (see reviews [1, 2]) have 
put graphene in the focus of current reserach in condensed matter physics and far beyond. The principal feature of 
graphene is the linear dispersion law u> = Vpk (vp — c/300 is a Fermi velocity, c is the speed of light in vacuum) 
for quasi-particle fermion excitations above graphene Fermi surface, which is valid for N — 4 species of fermion 
excitations in a graphene layer. This law is valid for the energies up to a few eV. The dynamics of quasi-particles and 
their interaction with electromagnetic field are therefore well described by the quasi-relativistic Dirac model [3] . The 
propagation of electromagnetic field in the presence of a graphene layer is governed by the effective action where the 
J^~j ■ fermionic excitations are integrated out. To the lowest order this action is determined by the standard polarization 
£■ — ' operator of 2 + 1 dimensional fermions. Such diagrams have been calculated in a number of papers, see [4-6]. One 
t-H , of the most spectacular predicitions of the Dirac model is the universal absorption law for light passing through 
suspended layers of graphene. This law was experimentally checked with a high precision [7]. There are also other 
experimental confirmations of the Dirac model, most important of which is the observation of the non-conventional 
H quantum Hall effect, [8]. 

t-H | In this paper we study the Casimir interaction at non-zero temperature of a graphene layer with a parallel metal, 
using the quasi-relativistic Dirac model for descriptions of quasiparticles in a graphene. In the framework of this 
model the problem was already addressed in [9] , where it was shown that at zero temperature the interaction between 
graphene and an ideal conductor is about 2.6% of the interaction between two ideal conductors separated by the 
same distance. Earlier, similar calculations were done [10-12] for the hydrodynamic model, which does not reproduce 
linear spectrum of low-energy excitations in graphene. More recently, thermal Casimir interaction between graphene 
layers was calculated in [13] in a van der Waals type approach, i.e. basing on the non-retarded Green's functions. 
In [14] in the framework of the quasi-relativistic Dirac model it was found that the Casimir interaction between two 
suspended graphene layers becomes strong at high temperature due to matter plasmonic response in a graphene layer. 
In another paper, Ref. [15], where graphene was modelled by a combination of Lorentz-type oscillators, practically no 
temperature dependence of the Casimir force was found. 

Temperature dependence of the Casimir interaction attracts much attention due to the unresolved problem of its 
asymptotics for metals at large separations and finite temperatures [16-21]. Two models of dielectric permittivity, 
the Drude model and the plasma one, lead to essentially different Casimir interaction due to their different infrared 
frequency dependence [22-25]. The free energy of two metals with parallel surfaces separated by a vacuum slit a, 
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for a > he/ (AitksT) is twice as much for the plasma model than for the Drudc one [22, 23, 26]. Various theoretical 
models were used to substantiate arguments in favor of the high temperature behavior characteristic for a Drude model 
[27-30]. Predictions of the models also differ at shorter separations where most experiments are being performed. 
However, a discussion of whether one of the models is excluded by existing experimental data is still under way [31, 32]. 
At the same time, we would like to emphasize that the main predictions of the paper are insensitive to particular 
model used for description of conductivity of metal (see discussion in Section IV. A). 

The results of the paper are outlined below in brief. In Section II we derive a polarization operator for 2+1 Dirac 
fermions at finite temperature and chemical potential. It allows us to describe in Sect. Ill optical properties of 
graphene in terms of the components of polarization operator using the modified Maxwell equations for classical 
electromagnetic field in the presence of a suspended graphene layer. Basing on this, we perform in Section IV a 
detailed study of the Lifshitz free energy of a graphene layer interacting with a parallel metal. In Sect. IV. A we 
derive the high temperature asymptotics of the free energy of the system, which appears to be the same as for two 
conducting surfaces described by the Drude model, i.e. surprisingly strong for a monoatomic graphene layer. In Sect. 
IV. B we investigate in detail the contribution of non-zero Matsubara terms, and discuss different scaling regimes of 
the free energy specific to considered graphene-metal system. In Conclusions we outline our main results. 

In what follows we adopt ^ = c = fcs = l. 



II. THE DIRAC MODEL 



The dynamics of quasi-particles in graphene is well described within a 2 + 1-dimensional Dirac model. All the 
most important properties of low-energy excitations are indeed readily incorporated in this model: the linearity of the 
spectrum at the energies below approximately leV, a very small mass gap (if any, but see discussion in [4-6, 33]), and 
the different "speed of light" inside the graphene layer vp ~ (300) -1 . Moreover, the interaction of the quasi-particles 
with cither quantum or classical external electromagnetic field is straightforwardly described within this model by 
a standard "covariantization" of the derivative with help of electromagnetic potential A^. Thus one deals with the 
model described by the following classical action (assuming graphene plane lying at x 3, — 0) 

S = ~\f d A xFl v + J d 3 x^ (1) 

with 

lp = (id -H- eA ) 7o + vfYj 1 (idx - eA^ + 7 2 (id 2 - eA 2 )} - m . 

Since there are N = 4 species of fermions in graphene, the gamma matrices are in fact 8x8, being a direct sum of 
four 2x2 representations (with two copies of each of the two inequivalent ones), Jq = — (7 1 ' 2 ) 2 = 1. The Maxwell 
action is normalized in such a way that 

47T 

6 = 4 ™ = 137' (2) 

The properties of the electromagnetic field (both classical, and quantum) in presence of a single graphene layer can 
be deduced from the partition function 



Z = j D[Atpi>]e iS . 



where S is given by (1). Depending on the problem under consideration we can also constrain the electromagnetic 
potential A with some additional conditions. In particular, investigating the Casimir interaction between graphene 
and parallel ideal conductor we will impose the conductor boundary conditions at x 3 = a 



i 0\x 3 =a 



M\x 3 =a = A 2 \ x 3 =a = dzM\x 3 =a = °- (3) 



thus making the partition function dependent on the distance between graphene and the metal, Z = Z(a). 
As the first step of consideration, one shall integrate out fermions thus obtaining 



/ 



At a somewhat formal level 

S eS (A) = -ilndet(^). (5) 
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Since the action (1) is quadratic in ip, the expression (4) is exact, though one has to give a precise meaning to the 
determinant (5). 

In the quadratic approximation the effective action S e s(A) discussed above reads 

1 f d3 P * f \~njl i 




S cS (A)=A~\ r A =2j -J^MpW{p)M P ), (6) 

The polarization operator Ily has been considered in a number of papers for the systems characterized by various sets 
of parameters. For the zero-temperature case it has been first calculated in [4], for non-zero T in [34], other cases were 
considered in, e.g., [6, 35]. Here we present the general result for Hij(po, p) depending on non- vanishing m, T . The 
final formulae (9), (13) coincide with the cases presented in the literature provided one considers appropriate limits 
of parameters. In particular, (A. 27) of [34] can be obtained from our result (13) for m = fi = and vf = 1, while 
(30) of [35] in the case of T = is equal to (13) in the limit po — > 0. The calculations are therefore quite standard 
though tedious. We shall be rather sketchy in their description. 

In Minkowski space the one-loop polarization operator can be expressed in momentum space as 

U mn ( Po , p) = ie 2 J tr (5(g , q)7 m 5( 9 o - Po, q - p)j n ) , (?) 



where the propagator of the quasiparticles in graphene reads 



a, v (Qo + Mho - vf$. - m 

(q + p + ie sgngor _ '"f^ ~ m 

Note that due to the quasi- relativistic nature of excitations in graphene S also depends on the Fermi velocity vp- 
Here q J = (q 1 ^ 2 ), $ = 7^1 +J 2 Q2- 

Due to the Lorenz and gauge invariance (transversality) the polarization tensor in an empty space can be completely 
defined by calculating a single scalar function. In the presence of medium characterized by the velocity u one needs 
two independent components for complete definition of the (parity-even) W J [36] 



,,2 'H 



nf;A(po,p)+p 2 Tli l B(p ,p) (9) 



jjji = ji _ P^f_ n ji = _ p>u l + 11^' v?u l _ 2 

p 2 ' u p 2 (pu) (pu) 2 

where 77 = diag(l, Vp, vp), and A, B are scalar functions. Note, that in this expression we also took into account the 
above mentioned quasi-relativistic nature of the excitation in graphene by introducing an appropriate dependence on 
Vp. Vectors with tilde are rescaled by multiplying the spatial components with vf, i.e., p 3 = v/lp 1 = (jpo,Vfp). For 
more details see, e.g., [37]. 

In the medium rest reference frame, u = (1,0,0), the scalar functions A, B can be expressed via the temporal 
component of the polarization operator and its trace in the following way (p 2 = p 2 , — p 2 , p 2 = p 2 , — v 2 F p 2 } 

A = ^|n o + n tr , b = ( P3+ 2 P3 n QQ + nj) . (10) 
p 2 ^p 2 vp / 

Clearly, we can express A, B through any pair of the components of n. We choose n tr = n™ and Hoo for convenience 
of further calculations. In what follows we sketch the derivation of the temporal component Hoo only, while the trace 
is calculated similarly. 

To introduce the temperature in (7) we perform the rotation to the Matsubara frequencies (see, e.g. [38]) 

/CO 
dq ^-2irT^2, q -> 2mT{k + 1/2), (11) 
k— — 00 

use the Feynman parametrization 

1 f 1 dx 

ab J (xa + (1 — x)b) 2 
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and subsequently change the variables in (7) in the spatial part of the loop-integration: q — > q + xp. Then we come 
to 



k= — oc 



d 2 q M 2 - ( 9o + m)(Po - 9o - m) 



) 2 [( go + /i-^o) 2 -e 2 
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with M5 2 = m 2 + uf,q 2 - a;(l - x)vf r p 2 and 

9 2 = m 2 + w 2 q 2 - x(l - x)(pl - v 2 F p 2 ). 
Summation over the fermion Matsubara frequencies can be made explicitly 

f " 2 = -77^ ( 6 sech 2 f 2 + - 2Ttanh ( *±±)) + (9 -> -9). (12) 

[(27riT(fc + 1/2) -fc) 2 -9 2 ] 2 169 3 T 2 V V 2T / V 2T y7 v y v ' 

The q-integration can also be performed directly by using that d x tanhx = sech 2 x. After these transformations only 
the integral over x remains, and we arrive at the following representation for n o and n tr 

n tr ,oo - I dx (f tr . 00 tanh - In U cosh + ( 0O -+ -Q )) (13) 



J F JO 



where 9o = \Jm 2 — x(l — x)(pf ) — v F p 2 ), b = p x — fi, and 

-2« 2 ,p 2 a;(l-x)-po(l-2.x)9o + 29 2 
/oo = ^ • (14) 

2m 2 v 2 F + 2x{l - x)v 2 F p 2 Po (l - 2v 2 F )(l - 2x) - 2(1 - v 2 F )O 
/tr _ 4T9^ 4T ' U5j 

We remind that N is the number of fermion species, N — 4 for graphene. Parity-odd contributions to the polarization 
tensor cancel out between different species, while the parity-even contributions add up. 

We note here that the expression (7) is power-counting divergent. An ultra-violet divergence indeed appears in the 
trace part II™. To remove this divergence we performed the Pauli-Villars subtraction at an infinite mass, which is a 
rather standard procedure. One can check that the m — > oo limit of (13) is indeed zero. 

III. REFLECTION COEFFICIENTS 

In our approximation the full action for the electromagnetic field is — j J d A xF 2 v + S e g, where the effective action 
(6) is confined to the surface of graphene, which we place at x 3 — 0. This action is singular and gives rise to an 
interaction proportional to S(x 3 ) and depending on the tangential momenta. Therefore, the Maxwell equations receive 
a singular contribution 

+ 5(x 3 )TT p A p = 0. (16) 

We extended II to a 4 x 4 matrix with n 3AI = IP 3 = 0. The equations (16) describe a free propagation of the 
electromagnetic field outside the surface x 3 = subject to the matching conditions 

"<4-/z|x 3 = + = A tl \ x 3 — _Q, 

(d 3 A fl )\ x3=+0 - (d 3 A lt )\ xa= _ = n»A v \ x3=0 (17) 

on that surface. So far, the only difference from the zero-temperature case [9] is in the form of U.J 1 . These matching 
conditions can also be rewritten in terms of the electric and magnetic fields E, H. 
Next, we introduce a TE mode 

E = (-p 2 ei + pie 2 )p *(x 3 ) (18) 
H = i(piei +p 2 e 2 )*'(x 3 ) + e 3 (p 2 + pj)^(x 3 ) (19) 

(20) 
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and a TM mode 

E = i( Piei + p 2 e 2 )$'(x 3 ) + e 3 (p? + P 2 2 )^(x 3 ) (21) 
H = (pad - Pl e 2 )p $(x 3 ) (22) 

(23) 

where ei, e 2 , e3 are unit vectors. An overall factor of exp(i(x°po + a^Pi + x 2 £> 2 )) has been omitted for brevity. 
To define the scattering data in the TE and TM sectors we take the potentials in the form 

(e^ + r TE e-^\ x* < + x * <0 

\tTEe tp3X , x 3 > y J \tT M e tp3X , x 3 > V ; 

The reflection and transmission coefficients are defined by the matching conditions. After some algebra we obtain 

A 2ip 3 
rTE = 2^4' = 2^4 (25) 



^2 



p 3 (A-^p 2 B) 2^ 

TTM — TTTTo 71 5 — 777": ITM 



2ipl-p 3 (A-v 2 F p 2 BY 2ipi- P3 (A-v 2 F p 2 B)- 

For the later use in the Lifshitz formula (27) for the Casimir energy, one can also rewrite the reflection coefficients in 
terms of the polarization tensor components 

P3n 00 _ pjiioo + p 2 n tr 

r ™~p 3 noo + 2z P 2 ' rTE " p|n o + p 2 (n tr -2 J p 3 )' [ ' 

The predictions of (26) for the visible light appear to be essentially the same as for the idealized (m = = 0) 
zero-temperature case [37] since the corresponding frequencies are much higher then any other scales in a realistic 
graphene sample. 

We stress that here we are working with free standing samples of graphene. Otherwise, the substrate will contribute 
to the reflection coefficients. 



IV. FREE ENERGY 



As early as in 1955 Lifshitz demonstrated [39] that the Casimir interaction between two parallel dielectric slabs can 
be expressed in a closed form if their dielectric permittivities are known at the imaginary frequencies. In a number of 
later works the original calculation was generalized and refined [40-43]. In particular, it was shown by Kats [40] that 
for any two given parallel plane interfaces separated by the distance a and described by their reflection coefficients 
r TE TM' r TE tm 01 tnc ^E and TM electromagnetic modes, the Lifshitz free energy density reads 

^ = T E /|^M(l-e-*'"rM)(l-e-*i-r«r«)]. (27) 

n— — oo 

where p\\ = y 'uj 2 % + p 2 , and w„ = 2imT are the Matsubara frequencies. The reflection coefficients here are assumed 
to be also taken at Euclidian momenta r — r(u n , p). 

Choosing as interacting interfaces a suspended graphene film and a parallel ideal metal, we shall substitute in the 
above the corresponding reflection coefficients. For the perfect conductor one has 

»™ = 1, ^te = -1- (28) 

The reflection coefficients for graphene at Euclidean momenta can be found after t he substi tution po = i2nnT = iui n 
into the formulas (13-15), (26). More precisely, in (26), one should replace p 3 by iy / w 2 + p 2 . 

In the following section V.A we consider the high temperature (large separation) asymptotics of the graphene-ideal 
metal free energy, while in subsequent section V.B we study the case m = jj, = analytically and discuss numerical 
results. 
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A. High temperature limit 

As discussed in the introduction, the high temperature limit is the most intriguing regime in the Casimir physics. 
In this limit at distances H = AiraT ^> 1, the zero frequency Matsubara terms determine the free energy since other 
Matsubara terms are exponentially suppressed at these separations. As we show below, at zero frequency ttm ~ 1 ; 
tte ~ OiNvp\p\ for small values of the wave vector |p|. Thus the coefficient ttm acquires the ideal metal value, 
and this is the principal feature of graphene at high temperatures. At the same time, the TE mode free energy for 
graphene- metal is suppressed at separations a > l/(4irT) by a factor aNvp / (aT) <C 1, and due to this suppression 
the difference in this regime in the results for the ideal metal and the Drude or the plasma models of a real metal is too 
small to be observed in possible experiments in graphene-metal system. For other distinctions of the graphene-metal 
high-temperature asymptotics see the next subsection. 

As already mentioned the main contribution in the sum in (27) at high aT is given by the zero frequency Matsubara 
term 

F Q = —j sds ln C 1 + ^T E e- 2aTs ) (1 - r TM e- 2aTs ) (29) 

where we also performed angular integration and re-scaled the momentum, s — \p\/T as compared with (27). Further 
expansion in powers of aT is determined by the s — > behavior of the integrand in (29). To obtain corresponding 
asymptotics we first note that the polarization tensor components at w n = behave as 



Hoc ^ -J3- In 2 cosh J - - tanh __£ + ->-„)] + (30) 



aNT 2 s 2 / , m 



m - 



+ 24m (2tanh^^ + — tanh^^^-m+Oi^-^)) +0(s 4 ) 



„ „ aNT 2 vis 2 ( , m + u , A .„„. 

n tr - n 00 ^ —^JL- (tanh -UL + (/1 _> _„) j + (s*) (31) 



From this we readily deduce that 



r TE (tanh + (m - -/*) j + 0(. 2 ), (32) 

, sv F /{aN) ,r>( 2 \ 

r ™s^o L m ( 2 cosh ^) -# tanh ^ + ( M ^- M ) +tA J ' ^ 

where we assumed that m and /i are of the order of the temperature T, which can be relevant for some particular 
graphene samples and/or experimental setups. Note that all formulas in this section are also valid for m, /j, — > 0, since 
s — > is essentially the same limit as T — > oo. 

As we can see, tte remains small for the relevant values of momenta. Taking the leading term in the expansion of 
the logarithm in (29) and using (32) we obtain the leading TE contribution at large separations corresponding to the 
zero frequency term 

*™ * " ( tanh "22^ + ^ "">) • (34) 

For the TM mode one cannot use the expansion of ln(l — r^M e ^ 2aTs ) as it is divergent near s = 0, see (33). Instead 
one has to use the representation 

*™ = g jf° sds (in (1 - e-™> ) + In (l - ^r^) ) • (35) 

The integral of the first term can be taken explicitly leading to 



-r-(O) _ k B T((3) _ _ x 

otm ~ 167ra 2 = Drudc l T ^°° — 2 id l T ^°°' ^ ' 
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which coincides exactly with the interaction of two real metals at high temperature described by the Drude model, 
or, equivalently, one half of interaction between two ideal metals at high temperature, J 7 ^. For readability, in (36) we 
temporarily restored the physical units. In the rest of (35) one can use an expansion 

4?M = -£gif ^(T^)'e-^ (37) 

As the main contribution to the asymptotics at H — » oo comes from vicinity of s = we can expand ttm — 1 in a 
Taylor series around this point, to obtain in the highest order 

T3 » i (r&Vg / / J_\\ ^ ^rg^i C(3)r& 

4tt^; (2aT)'+ 2 V \2aTjJ~ 4ir (2aT) 3 167ra 3 k j 

where C; = J °° s' +1 ( j-zpr; ) e?s and r^j = (J; r TM) s=0 which can be deduced from (33). In graphene, where 

vp/(aN) <C 1, this correction is already small for considered H = 4naT > 1, but in other Dirac systems with larger 
vp, the transition of J-otm (35) to asymptotic (36) would occur only at much higher distances aT ^> vp/{aN). We call 
attention of the reader that in any case this condition does not define the transition to the high-temperature regime 
of the whole free energy (27) which also includes the nonzero Matsubara terms. In the next subsection we investigate 
in detail their contribution and define through a scaling argument the transition condition for graphene-metal system 
as aT > alna- 1 /(2C(3)) ~ 0.015. 

Thus we can say, that from (34), (36) and (38), and considerations of the next subsection it follows that for 
separations H ^> 0.19 the main contribution to the interaction of a suspended graphene layer with ideal conductor 
comes from the TM mode leading contribution (36) and asymptotically constitutes just one half of the interaction 
between two ideal metals at large separations. This is induced by the specific dependence of the TE,M reflection 
coefficients on the momenta s = \p\/T for small s, see (32), (33). As mentioned earlier, such ttm,e behavior also makes 
the predictions for large separations insensitive to the differences between Drude and plasma models of conductivity of 
a real metal in possible experiments. A numerical analysis shows that the difference between graphene-ideal metal and 
graphene-gold free energies becomes small already at a > lOOnm at T = 300K (see Fig. 2) and rapidly decreases with 
higher separations. Moreover, in this limit the interaction of two suspended graphene samples will acquire exactly 
the same asymptotical value. 

As was shown in [9], in the zero-temperature limit the Casimir force between graphene and perfectly conducting 
metal is about 2.6% of the force between two perfectly conducting metals. As we see, in the opposite limit, H 3> 0.19, 
the ratio of graphene - ideal metal free energy to the ideal metal - ideal metal free energy is strongly enhanced. Such 
an enhancement takes place due to the non-perturbative structure of the Lifshitz result in the limit vp =0. 

B. Further analysis and numerical results 

Given the specific values of m and /i, the Casimir interaction can be evaluated by making use of the Lifshitz 
formula, reflection coefficients and the polarization operator. As our numerical analysis shows, increasing the mass 
gap diminishes the Casimir interaction (see also a corresponding discussion in [9]), while inclusion of a nonzero chemical 
potential fx enhances the interaction. However, the difference in the values of the free energy for m, fi ~ O.OlcV, which 
are reasonable bounds for these quantities in suspended graphene samples, and the m = fi = case is less than one 
percent. Therefore in what follows we turn our attention to the important case m = \i = and study it in detail. 

It is often desirable to have an accurate analytical approximation of the exact result at different separations. We 
present such an expression for the sum of nonzero Matsubara terms in this section. 

To obtain an appropriate analytical expression we first note that at separations H ^> vf one can put vf = in 
any nonzero Matsubara term. It is possible due to the exponential factor in the Lifshitz formula which effectively 
restrains the integration over impulse to apn < 1. In this case contribution of the type of vp(ap\\) 2 can be neglected 
compared to (ap n ) 2 = {2imaT) 2 due to the smallness of the parameter vf- However, our numerical analysis shows, 
that this approximation works at smaller distances as well, see below discussion around (46). 

Secondly, in the finite temperature sum of nonzero Matsubara terms in the Lifshitz formula one can use the 
reflection coefficients taken at zero (!) temperature. The corrections due to finite temperature are suppressed for 
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FIG. 1: Ratio pi of the free energy to the high-temperature asymptotics Jq T m, see Eq. (36). Both graphs are calculated for 
T = 300K. In graphene the values m = fi = were used. 



nonzero Matsubara terms, so we neglect them in the leading approximation. Contrary to this, in zero frequency 
Matsubara terms one has to take into account the whole structure of the finite temperature polarization operator 
entering ttb,M- 

Under two mentioned above approximations and the condition m = /J, = the reflection coefficients of a single 
graphene layer at zero temperature have the form: 



T _ naN^Jpl + p 2 TraN^Jp 2 +p 2 

r TM = / / = — / (39) 

iraNJp 2 , + p 2 + 8^jp 2 + v 2 F p 2 iraN ^jp 2 + p 2 + 8 \p 1 

r T=o_ ™N^pI + vIpI ^bol (m 

' TE — I = 1 — / Y^) 

naN ^Jpl + v 2 F p\ + 8^Jp 2 + pjj iraN\p \ + 8 ^Jp 2 + pjj 

Due to smallness of the reflection coefficients (both being of the order of a) we can take just the first term in the 
expansion of the logarithm in the Lifshitz formula as another reasonable approximation for nonzero Matsubara terms. 
Note, however, that expansion of the reflection coefficients themselves (at least in TM mode) is not legitimate, as will 
become evident below. 

The sum of nonzero Matsubara terms in this approximation in the TM case with rl 1 ^ = t"tm° (^9) an d r TM = ^ 
equals to 

AJ"tm = -^-V/ dsi — — ^— — --^exp(-2asi) = 



rp +00 



8na 



exp(-Hn)(l - gn + Hn) + (gn) 2 exp(gn)Ei(gn + Hn), (41) 



here Si = yj ui 2 + p| and g = 32Ta/(aN), E\ stands for the standard exponential integral function. It is convenient 
to reexpress the result (41) in an integral form. For this purpose one has to differentiate (41) over H , assuming H as 
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FIG. 2: Ratio p2 of the free energy for a graphene - metal system with fi = m — to the ideal metal - ideal metal free energy 
at T = 300K. 



an independent parameter for the moment, calculate the sum over n and then integrate back over H (the integration 
constant is fixed as zero at H — > oo). Thus one obtains 

ATxm = -™ f°° eft -PW(^) + 1 )^ 2 (42) 
Jff (exp(t) - 1) (8iJ + naNt) 

At large separations one can further approximate 8H + iraNt ks 8H (i.e. expand in powers of a) and get 

aN ( , . , Hexpt-H) iJ 2 exp(-F) \ _ 2n 

lzona^ \ (1 — exp(— H )) 2(1 — exp(— rl)y J 

The formula (43) provides an accurate description of (42) for H > 1 where AJ-tm constitutes less than 6% of the 
whole answer. At T — 300K this condition corresponds to distances greater than 600 nm. For smaller separations it 
is convenient to use (42) directly. 

The TE part of the nonzero Matsubara terms of the Lifshitz formula with the coefficients r!p E = Tte° from (40) 

(2) 

and r T ' E = —1 gives the following contribution in the leading order in a 

T 2 TraN±^ f +0 ° , , . T 2 TraN exp(-if) 

AF TE ~ — d Sl exp -2a Sl = ttt— j, U , ' ■ 44 

8 J Un 16a (1 -exp(-H)) 2 

Thus, the complete result for the sum of nonzero Matsubara TM and TE terms in the approximation described 
above is given by (42) and (44) 



AT = AJtm + AT TE - (45) 

Consequently, the leading vp = contribution to the free energy is the sum of (36) and (45): + AT. It can 

be used for the comparison of the theory and experiment almost at all separations. As our numerical analysis shows, 
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a difference between the exact Lifshitz result T (27) (with vp = 1/300, m = \x = 0) and considered approximation 

•^otm + i s l css than 1% for all separations between 1 and 2000 nm at the temperature T — 300K, which 
approximately corresponds to 3.3 > H > vp/2 — 1/600. 

In our opinion, the distances corresponding to H > vf (i.e. above ~ 2 nm for T = 300K) are most relevant for 
a possible experiment, and the analytical expressions (42) and (44) are obtained to describe this particular regime. 
However, numerical studies reveal that £±Ttm,e are smooth in the vicinity of H = point, and collate with T = 
result obtained in [9]. Indeed, from (42) and (44) one gets in this limit 

^^r-T^M^^N*))-^. (46) 

where the non analyticity in a comes from the TM mode, and explains the precautions we took avoiding to expand 
the reflection coefficient in (39). Since for T — > the zero- frequency term (29) vanishes, eq. (46) represents the 
free energy in the leading order in a in the approximation when vf — 0. The difference between (46) and the exact 
Lifshitz result at T = with vf = 1/300 is less than 1%. We also note that at T = the perturbative expansion in 
a is meaningful for vp ^ only, see [9]. 

On Fig.l we plot the ratio p\ of the Lifshitz free energy T (27) divided by the high temperature asymptotics J^tm 
(36) at T = 300K for two systems: graphene- ideal metal and ideal metal-ideal metal. From Fig.l it is evident that 
the pace of approaching the asymptotic values by the free energy in the case of a graphene - ideal metal system is 
much higher than in the ideal metal-ideal metal one. It is explained by the particular behavior of the components of 
the polarization operator at nonzero frequencies 

n oK) ? ^n + o(p 4 ), n tr (o ^ aT~c n + o( P 2 ), n >i, 

|p|-s-0 i |p|->o 

which induce additional 0(a) suppression of the corresponding contributions to the free energy. Here, c„, c n are 
factors of the order of unity. Because of this suppression the zero frequency TM Matsubara term dominates in the 
temperature dependence of the Lifshitz free energy for a graphene-metal system at much shorter separations than in 
the case of a metal-metal system. 

It is instructive to define separate regimes on the basis of the scaling behavior of the free energy, rather than on 
validity of various approximations. As a function of distance, T ~ a~ 3+s , where 5 = at T = and 5 = 1 at T — > oo. 
The scaling parameter 6 varies slowly, and can thus be defined through a logarithmic derivative 

S = a^-[\na 3 \T\] . (47) 
Since we are interested in rough estimates, it is enough to approximate J 7 by a rational function 

x- U + aV , Aa \ 

^-le^r- (48) 

where U can be defined by the T = value (46) and V is calculated through the high temperature asymptotics, (36). 
Namely, 



In 1 



1 



aNn J 2 



V = TC(3) . (49) 



Consequently, (47) gives 



" T -mh-y < 50) 



We may consider the "zero temperature" regime as the one where 5 is close to its T — value, let us say 5 < 0.2, giving 
aT < 0.004 or a < 30nm at T = 300K. The high temperature regime, on the other hand, may be defined as S > 0.8 
corresponding to aT > 0.06, or a > 460nm at T = 300^. Numerically, this bound is very close to one in the original 
condition, H ^> 1, which is universal for any Casimir system. However, in the graphene-metal interaction, as we 
already discussed, the high-temperature asymptotics is already saturated at H ~ 1, see Fig. 1. This can also be seen 
by considering an 'intermediate' scaling. At distances corresponding to S = 1/2, where H ss 27ralna~ 1 /((3) w 0.19 is 
much smaller than unity, the high temperature asymptotic value already becomes larger than zero-T contribution, i.e. 
aV > U, and thus this point can be considered (rather formally) as a crossover between zero and high-temperature 
regions. 
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It is interesting to note that the parameter U which governs transitions between different regimes is defined by a 
and does not depend on vp. One should not however overestimate this fact. All our expressions are valid for small 
vp only. In such a case the free energy behaves, very roughly, as a sum of the zero temperature contribution and 
the zeroth Matsubara term. Therefore, the crossover between different regimes is defined by the ratio of these two 
contributions. 

On Fig. 2 we plot the ratio p2 of the free energy of a graphene-metal system to the ideal metal-ideal metal free 
energy at separations below 300 nm. The temperature is chosen as T = 300K, which means that H — A-KaT varies on 
Fig. 2 from zero to 0.49. We used a plasma model for gold e(icu) = 1 + uj p /uj 2 with a plasma frequency uj p — 9.0eV in 
numerical calculations of the free energy for gold semispace and a parallel graphene layer, gold reflection coefficients 
are standard Fresnel coefficients for TM, TE modes with a dielectric permittivity e(w), see [44]. The numerical 
result for the free energy of gold-graphene system divided by an ideal metal - ideal metal free energy is shown by 
the dashed blue curve on Fig. 2. Note that the behavior of the gold-graphene free energy at short separations is 
one power of a different from the graphene-ideal metal free energy. This is the usual change in power law when 
one approaches short separations and the transition from the retarded regime to a non-retarded one takes place at 
separations characterized by a material wavelength X p = 2tt/cj p . Thus, the result for a real metal interacting with 
graphene is essentially different from the ideal metal - graphene result at short separations only. The red curve which 
corresponds to the zero Matsubara term only is almost perfectly linear. This reflects the point that, as mentioned 
above, the expansion of J-qtm around the high temperature asymptotics ^qtm ^ s governed by a small parameter 
vp/(aN), while the term J-qte is suppressed by a factor aNvp. 

Physical units can be restored in our formulas following the simple rule: in the result expressed as an appropriate 
power of aT divided by a 3 , the former product is to be replaced by fc^aT '/ '(he), while 1/a 3 is to be substituted by 
he /a 3 . 

V. CONCLUSIONS 

We have calculated finite temperature Casimir interaction between graphene and a parallel conducting plane in the 
framework of the Dirac model of quasiparticles in graphene. From the theoretical point of view, it is interesting to 
note, that we start with a fully consistent quantum field theory model at finite temperature, and that the polarization 
tensor (conductivity of the graphene surface) is thus temperature dependent. 

At high temperature (large separations) the free energy of interaction of suspended graphene sample and a parallel 
metal asymptotically behaves as the one for two conducting surfaces described by the Drude model, which yields a 
very strong Casimir interaction, somewhat surprising for a one-atom thick system. This feature provides an excellent 
opportunity for the experimental studies of the temperature Casimir effect in this system at room temperature and 
large separations as well as for medium ones. The energy at large separations is essentially insensitive to the model 
which was actually used for description of the metal's conductivity. Moreover, the interaction between two graphene 
samples would have exactly the same hight temperature asymptotics. Such behavior is induced by the specific static 
limit (w — 0) of the graphene reflection coefficients, rpu — I^te — given the small Fermi velocity in graphene, 
Vf <C 1. In this limit they formally coincide with the reflection coefficients for the Drude model of metals. 

We studied the Lifshitz free energy for graphene-metal system in detail and obtained results that can be readily 
used for the comparison of the theory and experiment at room temperature and various separations. For separations 
corresponding to H = AiraT > vp, we used the limit vp = in the Lifshitz formula to evaluate the sum of nonzero 
Matsubara terms. The speed of quasiparticles in graphene is much slower than the speed of photons, this is the 
physical reason for the accuracy of the vp = limit. In this limit perturbative results in the coupling constant a 
diverge at T = 0. Perturbation theory in a does not give reliable results for graphene systems where vf is small. 

Different physical regimes were separated by investigating the dependence of the free energy on the distance, a. 
It appeared to be close to the typical zero temperature one, a -3 , for ksTa/ (he) < 0.004 and approaching the high 
energy asymptotics a~ 2 for fcgTa/ (he) > 0.06. The formal crossover between different regimes can be chosen as 
kBTa/(hc) ps 0.015. It is interesting to note that the metal-graphene system arrives at the high-temperature regime 
much faster than the ideal metal-ideal metal system due to the 0(a) suppression of the TE mode and all nonzero 
Matsubara terms. 

We have to stress that we have considered a suspended graphene sample. In such a sample the values of m and fx are 
small and, as our numerical study shows, their influence on the Casimir interaction is negligible. This approximation 
will not be applicable to graphene on a substrate, where one should also take into account the influence of impurities 
and the contribution of substrate itself to the reflection coefficients. Another point of special attention which should 
be taken into account when comparing the theory with experiments is the fact, that exact results in complicated 
geometries [43, 45-48] can be essentially different from the approximations based on the Lifshitz result for two parallel 
plates. In any case, the state-of-the-art experimental techniques resolve the Casimir force with total error of fractions 
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of percent. This permits us believe that a graphene based experiment being not an easy task is still perfectly feasible. 

Finally a comparison with some previous works on the subject is in order. As mentioned before, the leading 
asymptotics for the free energy of a graphene-metal system at large separations (36) coincides with one of a graphene- 
graphene system. Thus our findings support the results of Ref. [14] where the thermal van der Waals interaction in the 
later system was studied. However, our prediction for characteristic distance separating zero and high-temperature 
regimes (as discussed below Eq. (50)) differs from that of Ref. [14]. In our case it does not depend on the Fermi 
velocity vf but does depend on the fine structure constant a, while in [14] it is proportional to vp with no dependence 
on a. There is no direct contradiction, however, since different physical systems have been studied in our works. 
We disagree with the paper [15], which found practically no temperature dependence of the Casimir interaction of 
graphene (and which therefore also contradicts [14]). Still, the estimate for the zero-temperature case given in [15] 
coincides with our previous calculations [9]. Finally, the separations considered in [13] are too small to allow for a 
comparison with our results. 
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